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Abstract 



In |[T1 an integrable quantum model was introduced and a class of its cyclic representa- 
^ ■ tions was proven to define lattice regularizations of the Sine-Gordon model. Here, we 

^ \ analyze general cyclic representations of this integrable quantum model by extend- 

\0 \ ing the spectrum construction introduced in []2l in the framework of the Separation of 

^ , Variables (SOV) of Sklyanin. We show that as in [jT] also for general representations, 

O ■ the transfer matrix spectrum (eigenvalues and eigenstates) is completely character- 

ized in terms of polynomial solutions of an associated functional Baxter equation. 
Moreover, we prove that the method here developed has two fundamental built-in fea- 
^ \ tures: i) the completeness of the set of the transfer matrix eigenstates constructed from 

\ the solutions of the associated Bethe ansatz equations, ii) the existence and complete 

characterization of the Baxter Q-operator. 
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1. Introduction 

In [HII an integrable quantum model was introduceci^ and a class of its cyclic representations 
was proven to define lattice regularizations of the Sine-Gordon models in a given sector of 
the quantum theory. Main results of that paper are the simplicity and the completeness of the 
characterization of the transfer matrix spectrum. 

Here, some comments are probably in order to better point out the relevance of these results. 

Integrability: It is of common use to call integrable a quantum model which admits a one 
parameter family of commuting conserved charges. This is the case in the framework of the 
quantum inverse scattering method [TOl |3l where the transfer matrix defines the family of 
commuting conserved charges. Anyhow, it is worth remarking that a more precise definition 
of quantum integrability requires that this family of commuting charges is a complete set of 
observables; i.e. its spectrum has to be simple (non-degenerate). This last statement has to be 
independently proven as it was done in [[Tl for the class of cyclic representations there analyzed. 

The first fundamental task to solve for a given integrable quantum model is the exact solution of 
its spectral problem, i.e. the determination of the eigenvalues and the simultaneous eigenstates 
of the complete set of its commuting conserved charges. 

Completeness of Bethe Ansatz: In the framework of the Bethe ansatz several methods has been 
introduced to analyze this spectral problem which lead to the characterization of the eigenvalues 
(and for some of those to the eigenstates) in terms of solutions to an associated system of Bethe 
ansatz equations. Hereafter, a rigorous proof of the completeness of the spectrum characteriza- 
tion was a fundamental and complicate long standing goa|^in the literature. In yj this result has 
been reached showing in particular that the eigenstates of the transfer matrix constructed from 
the solutions of an associated system of Bethe ansatz equations form a basis of the representa- 
tion. Note that early there were only a few examples of integrable quantum models where the 
completeness has been proven, including the XXX Heisenberg model; see [ | 14|1 and references 
therein. 

There is a double motivation for the present article. On the one hand, we are interested in 
showing that the strong results of simplicity and completeness of the transfer matrix spectrum 
can be proven for general cyclic representations of the SG model. On the other hand, we want to 
use these representations to define a method to characterize the spectrum which is based only on 
the Separation of Variables (SOV) of Sklyanin [[TSl [T6l [T71 . Then, these general representations 

'From here on we will refer to it as SG model. 

-This lattice regularization goes back to ||3] |4] |5] and it is related to formulations which have more recently 
been studied in IHEllll. 

^It was addressed also by numerical methods like the one introduced in lfm[T2ll . see also |fT3]| for an application 
of it. 
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are of special interest as they define concrete examples of cyclic representations of quantum 
integrable models to which the standard construction^ [|24l|- ll28l of the Baxter Q-operator by 
using cyclic dilogarithm functions does not apply in general. 



1.1 Methodological aim 

In the framework of quantum integrability, there are several methods to analyze the spectral 
problem as the coordinate Bethe ansatz [|29l |30l |3T1, the TQ method [[30l . the algebraic Bethe 
ansatz (ABA) [|9l[T0l[3l, the analytic Bethe ansatz ||32ll . However, they suffer in general from one 
or more of the following problems: i) Reduced applicability; i.e. there exist important examples 
of quantum integrable models to which these methods do not apply, ii) Analysis reduced only 
to the set of eigenvalues; i.e. they do not allow for the construction of the eigenstates. iii) 
Lack of completeness proof, i.e. the completeness of the spectrum description is not assured by 
the method while there are known cases for which these methods lead to incomplete spectrum 
description [|33l . 

The separation of variables (SOV) method of Sklyanin is a more promising approach. Indeed, 
the SOV method resolves problems like the reduced applicability of other methods; in particular, 
it works for integrable quantum models to which ABA does not apply. Moreover, it leads to 
both the eigenvalues and the eigenstates of the transfer matrix with a construction which has as 
built-in feature the completeness of the spectrum description. 

In the case of cyclic representations [34J of integrable quantum models, however, it is worth 
pointing out that the SOV method leads to the characterization of the transfer matrix spectrum 
in terms of solutions of an associated finite system^ of Baxter-like equations. Thus, SOV method 
does not lead directly to the standard spectrum characterization expected in the framework of 
Bethe ansatz approaches. Moreover, such SOV characterization of the spectrum is not the most 
efficient in view of the analysis of the continuum limit. In order to solve these problems it is 
important to prove that the SOV characterization can be reformulated in terms of an associated 
functional Baxter equation. One possibility is to add to the SOV characterization the Q-operator 
approach. In particular, this can be done when the corresponding Baxter equation is proven to 
be compatible with the finite system of Baxter-like equations of the SOV characterization. In 
fact, this was the strategy followed in [[T] for the special cyclic representations there analyzed. 

The main methodological aim of the present article is to present an approach which allows to 
reduce the solution of the spectral problem to the classification of the solutions of an associated 

'^Here, we refer to the Baxter's construction by gauge transformations which leave unchanged the transfer 
matrix while make triangular the action of the Lax matrices on the Q-operator. Let us recall that there are also 
others constructions of the Q-operator; interesting examples are presented in llT8l - [l23l . 

^The number of equations in the system is finite and related to the dimension of the cyclic representation. 
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Baxter equation in a fixed class of functions. An approach where the SOV characterization of 
spectrum is the starting point and the standard construction by gauge transformations§ of the 
Q-operator is completely bypassed. Let us recall that such an approach has been first developed 
in [|2l for the special class of cyclic representations studied also in [IJ. However, in [[H we 
have used some shortcuts to define the coefficients of the associated functional Baxter equa- 
tions which hold only for such a special class of cyclic representations. Here, we present this 
approach in the framework of the general cyclic representations also to infer some universal 
features of the method for applications to others integrable quantum models. 

1.2 Organization of the paper 

In Section 2, it is recalled the definition of the SG quantum model in the framework of the 
quantum inverse scattering method. In Section 3, the construction of the SOV representations 
corresponding to general cyclic representations of the SG model is implemented. There, the 
result derived in [[T] is extended and completed with the characterization of the coefficients of the 
SOV representations. Section 4 is the core of the article; there it is explained the method which 
allows to reformulate the SOV characterization of the spectrum in terms of an associated Baxter 
functional equation. Moreover, there it is proven the existence of all the polynomial solutions 
of this equation required to completely characterize the transfer matrix spectrum (eigenvalues 
and eigenstates) of the SG model. Appendix A and B contains some important properties. In 
Appendix C, the comparison between our SG model and the r'^^^-model is made, pointing out 
the difference in the spectrum of the two models. In particular, it is shown as for even chains 
the transfer matrices are not similar for general representations. 

Acknowledgments. I would like to thank J. Teschner for stimulating discussions on related subjects and 
for the opportunity to collaborate with him at DESY during the last three years under support of his grant 
Marie Curie Excellence Grant MEXTCT-2006-042695. I would like to thank also B. McCoy and J.-M. 
Maillet for the interest shown in this work. I gratefully acknowledge my current support from National 
Science Foundation, grants PHY-0969739. 



^Note that this method can be applied only when the existence of some model dependent quantum dilogarithm 
functions Il35l - ll45l is proven. 



2. Lattice SG quantum model 
2.1 Definitions and first properties 

The Lax operator of the SG model, as introduced in fll], reads: 



where A„ = A/^„ for any n G {1, N} with ^„ and parameters of the representation. For 
any n E {1, N} the couple of operators (u„,v„) define a Weyl algebra W„: 

u„v^ = g''"'"v^u„, where q = e-^'^\ (2.2) 

We will restrict our attention to the case in which g is a p-root of unity: 

/J^ = p = 2/ + l,p' = 2/' and /,/'gZ>° ^ = I. (2.3) 
P 

In this case each Weyl algebra Wn admits a finite-dimensional representation of dimension p. 
Let us denote with: 

I z ) = I zi, . . . , ) with G §p = = 0, . . . , 2/} and i G {1, N}, (2.4) 

the generic state of the basis constructed by the tensor product of the eigenstates | ) of each 
operator Un, then in this basis the Weyl algebra representation reads: 



(2.5) 



V„ I Zl, . . . , Zn) — t>„| 2^1, . . . , qZni • • • , Zf^) , 
Un I 2^1, . . . , Z^) = UnZn\zi,..., Z^) . 

where n„ and Vn are other parameters of the representation. 

The monodromy matrix of the model is defined in terms of the Lax operators by: 

and it satisfies the quadratic relations: 

R{X/fi) (M(A) ® 1) (1 ® M(/i)) = (1 ® M(^)) (M(A) ® l)R{X/fi) , (2.7) 
w.r.t. the six-vertex i?-matrix: 



f qX-q-^X-^ \ 
A-A-i q-q-^ 
q~q~^ A-A-i 
V qX- q-'X-' J 



(2.8) 
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Then the elements of M(A) generate a representation 7?.n of the so-called Yang-Baxter algebra 
characterized by the 4N parameters k = . . . , kn), ^ = (^i, . . . , ^n), u = {ui, . . ., mn) and 
V = (f 1, . . . , wn). In particular, the commutation relations (|2.7I) lead to the mutual commutativ- 
ity of the elements of the one parameter family of operators: 

T(A) = trc2M(A), (2.9) 

known as transfer matrix. 

In |[T1 the spectral problem of this transfer matrix has been solved for untwisted representations: 
kI e m, e M, < = 1, = 1, Vn G {1, N}. (2.10) 
In the present paper we will extend the analysis considering general representations: 

KleR, e M, (m„)*m„ = 1, (t;„)*t;„ = l, Vn e {l,...,N}, (2.11) 
which we call twisted for u"^ ^ 1 and v^f ^ 1. 

Let us comment that the above restrictions on the parameters of the representation are imposed 
to get the self-adjointness of the transfer matrix: 

Lemma 1. If the parameters of the representation satisfy the constrains (|2J7]) and 

e = —{i^nin)l {i^nCi) Uniform along the chain, (2.12) 

then the generators of the Yang-Baxter algebra satisfy the following transformations under 
Hermitian conjugation: 

, , ( Bt(A) \ ^ ( D(A.) C(.A.) \ 

^ ' \ Ct(A) Dt(A) ) \ B(£A*) A(A*) ) 

which, in particular, imply the self-adjointness of the transfer matrix T{X)for real A. 

Proof. It is simple to observe that the Lax operator of the SG model satisfies the equation (12.131 ). 
which can be also written as: 

(Lf(A))t = aiLf(5A*)ai ^ (M(A))t = aiM(5A*)ai, (2.14) 
that is (12.131) holds, when one takes into account that A(A) and D(A) are even in A. □ 
In the case of a lattice with N even quantum sites, we can introduce the operator: 

N 

e = nvi-^^"", (2.15) 

71=1 

which plays the role of a grading operator in the Yang-Baxter algebra: 
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Lemma 2. (Proposition 6 of ^) Q commutes with the transfer matrix and satisfies the fol- 
lowing commutation relations with the entries of the monodromy matrix: 

eC(A) = gC(A)0, [A(A),e] = 0, (2.16) 
B(A)e = geB(A), [D(A),e] = 0. (2.17) 



Moreover, the 9-charge allows to express the asymptotics of the transfer matrix as: 



N 



lim A±^T(A) = I n 1 (0 + • (2-18) 



Let us denote with Sj the set of the eigenvalue functions t(A) of the transfer matrix T(A). By 
the definitions (|2.1I) and (|2.9I) . Sj is a subset of ]R[A^, A"^]n/2' where we are using the notations: 

N = N + cn - 1, and eN = for N odd and 1 for N even, (2.19) 

and M[a;, x'^jn denotes the linear space in the field M of the real: 

f{x)e^x,x-^]M ^ (/(a;))* = /(x*) Va: G C, (2.20) 

Laurent polynomials of degree M in the variable x. 

Note that in the case of N even, the 9-charge naturally induces the grading Ej = UfcLo ^t'^' 
where: 

S^''^ = |t(A) G St : ^^^fiin ^ A±^t(A) = ^J] j {q^d + {q^d)-^) | . (2.21) 

This simply follows by the asymptotics of T(A) and by its commutativity with 0. In particular, 
any t{\) G Sj'' is a T-eigenvalue corresponding to simultaneous eigenstates of T(A) and 6 
with 9-eigenvalue q^9. 



3. Cyclic SOV representations 

The separation of variables (SOV) of Sklyanin [fTSl - lfTTl is a method to solve the spectral prob- 
lem for T(A) which is based on the observation that such a problem simplifies considerably if 
one works in a SOV representation: A representation where the commutative family of opera- 
tors B(A) is diagonal. This is due to the simple form assumed by the operator families A(A) and 
D(A) in such a representation. 
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3.1 Cyclic SOV representations: Generality 

In subsection l3.3l we will show that the commuting family of operators^ Bn(A) is diagonalizable 
with simple spectrum for almost all the values of the 4N parameters k, ^, u and v of our model. 
The corresponding SOV representations of the SG model are to a large extend determined by 
the Yang-Baxter algebra (12.71 ); here we recall the general form of these representations. 

Let ( I be the generic element of a basis of eigenvectors of Bn(A): 



6, (A) ^ n 



n=l a=l 



(3.1) 



where [N] = N — and 



(3.2) 



Here, the simplicity of the spectrum of Bn(A) is equivalent to the requirement ^ Q for any 
a 7^ 6 G {1, . . . , [N]}. The remaining generators of the Yang-Baxter algebra read: 



An(A) = e^\{\) 
Dn(A) = e^\{\) 



' N \ ' 



N 



At- - !^TH 



+ 2^ [[—7- -l—'^^iVa) T+ 



a=l b^a 



Va/Vb - Vb/Va 



.Vd " ^ 
where T^ are the operators defined by 

{r]i,...,T]^\T^ = {r]i,..., q^^ria, . . . , r/N | • 
While, Cn(A) is uniquel>§ defined by the quantum determinant relation: 

detqM(A) = A(A)D(g-U) - B{X)C{q'^X), 
where detqM(A) is a central elemenl§of the Yang-Baxter algebra (12.71 ) which reads: 



(3.3) 
(3.4) 



(3.5) 



N 



detqM(A) = Yl «^n(-^//"n,+ " Ain,+/A) (A//i„ _ - /i„, -/A), 



(3.6) 



n=l 



^From here, we will use the index N when it will be need to point out that we are referring to the chain with N 
sites and we will omit it otherwise. 

'^Note that the operator Bn (A) is invertible except for A which coincides with a zero of Bn, so in general Cn (A) 
is defined by (4.5) just inverting Bn(A). This is enough to fix in an unique way the operator Cn being it a Laurent 
polynomial of degree [N] in A. 

''The centraUty of the quantum determinant in the Yang-Baxter algebra was first discovered in P6) : see also 
BTI for an historical note. 
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where = ±m^^g^/^^„. Remark that the quantum determinant does not depend from the 
parameters m„ and Vn of the representation and so it is the same for twisted and untwisted 
representations. 

The expressions (13.31 ) and (13.41) contain complex- valued coefficients r]f^, rj^, aN(?7r) and d^{rir). 
Note that Za = r]^ and Zq = i]^ are fixed by the asymptotics relations (|B.9I) while the coeffi- 
cients a^{r]r) and d^{r]r) are restricted by the quantum determinant condition: 

det^M{r]r) = M^r)d^{q-\) , Vr = 1, . . . , [N] . (3.7) 

In a SOV representation, some freedom is left in the choice of ajsiiVr) and djq{r]r) which can be 
parametrized by what we call a gauge transformation: 

which just amounts in a renormalization in the states of the B-eigenbasis: 

N 
r=l 



3.2 Central elements: Average values 

Following ll34l . let us define the average value O of the elements of the monodromy matrix 
M(A) as 

p 

0{A) = l[0{q''X), A = A^ (3.10) 

k=l 

where can be An, Bn, Cn or Dn and we have to remark that the commutativity of each family 
of operators An (A), Bn(A), Cn(A) and Dn(A) implies that the corresponding average values are 
functions of A. So that ^n(A), 2^n(A) are even Laurent polynomials of degree N while -Bn(A), 
Cn(A) are odd Laurent polynomials of degree [N] in A. 

Proposition 1. 



a) The average values ^n(A), i3N(A), C^{h), 'Dn(A) of the monodromy matrix elements are 
central elements which satisfy the following relations: 

(^n(A))* = Pn(A*), (Sn(A))* = Cn(£A*), (3.11) 



under complex conjugation. 
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b) Let A^(A) be the 2x2 matrix with elements the average values of the elements of the 
monodromy matrix M(A) , then it holds: 

M^{k) = C^iA) C^^,{A) . . . C,{A) . (3.12) 

where CniA) is the 2x2 matrix with elements the average values of the elements of the Lax 
matrix L^'^{X). 



A similar statement was first proven in [1341 . 

Proof of a). Centrality of -Bn(A) trivially follows from the fact that Bn(A) is diagonal in the 
SOV representation and from (|3.1I ) it is easily found: 



N 



i3N(A) = z^^- n 



(A/Za-ZJA), _^ (3.13) 

n=l ' a=l = • 



The requirement of cyclicity of the SOV representation reads: 

(T:r = (T:)^ = l Vae{l,...,N}, (3.14) 

so that A-^{Zr) and 'D-^{Zr) are centrals and related to the coefficients a^^{q'^rir) and d^{q^r]r) 
by 

p p 
^N(^r) = n"N(gV), T^^{Zr) = JJ^NlgV), Vr G { 1 , . . . , [N] } . (3.15) 

k=l k=l 

Note that ^n(A)A^ and Pn(A)A^ are polynomials in A^ of degree N. So the centrality of 
^n(A) and T>f^{A) follows from the simplicity of the BN-spectrum and the centrality in the 
special values (13.151) to which we have to add for even N the centrality of the leading asymptotic 
terms of ^n(A) and I'm (A) as discussed in appendix |Bl Finally, the Hermitian conjugation 
properties of the elements of the monodromy matrix: 

(AN(A))t = Dn(A*), (BN(A))t = Cn(£A*), (3.16) 

imply (|3.1 II) and in particular the centrality of Cn (A) . □ 

Proof of b). Under the assumption that B(A) is diagonalizable and with simple spectrum in 
the entire chain as well as in each subchain, the point b) follows inductively by using the next 
Lemma. □ 

Lemma 3. (Proposition 3 of lH]]) The following recursive equations on the average values 
hold: 

B^iA) = ^M(A)i3N-M(A) + SM(A)PN-M(A), (3.17) 
Cn(A) = Pm(A)Cn-m(A)+Cm(A)^n-m(A), (3.18) 
^n(A) = ^m(A)^n-m(A) + Sm(A)Cn-m(A), (3.19) 
Pn(A) = Pm(A)I)n-m(A)+Cm(A)Sn-m(A), (3.20) 
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where on the l.h.s. there are average values of the monodromy matrix elements on the com- 
plete chain with N-sites while on the r.h.s. there are those of the monodromy matrices on the 
subchains i and 2 with (N — M)-sites and yi-sites, respectively. 

Remark 1. Being the average values central elements of the representation they are unchanged 
by similarity transformations and they therefore represent parameters of the representation. 
Moreover, the gauge transformations clearly leave Zr, Zi\, Zo, ^^(^r-) and 'D^(Zr) unchanged. 
Therefore these last numbers characterize gauge-invariant dates of the SOV representations. 



3.2. 1 Calculation of average values 

It is a simple exercise to show that the average values of the elements of the Lax matrices i^^'^(A) 
are explicitly given by 

where we have used the notations Kn = k^, X„ = in^Un = u^^ ^'^d = ^n- This formula 
together with equality (13.121) allows to uniquely define the average values of the monodromy 
matrix elements and to state: 

Lemma 4. Oj^(A)/ nn=i ^nfor O = A,B,C,V are Laurent polynomials of maximal degree 
1 in each of the parameters Kn, Xn, Un, Vn. 

Note that the previous lemma also implies that the gauge-invariant dates of the SOV represen- 
tations, Zr, Zj^, Zd, Aj^{Zr) and Vj^{Zr) (up to permutations of r = 1, . . . , [N]), are uniquely 
defined in terms of the parameters of the representation Kn, Xn,Un,Vn. 



3.2.2 Choice of the gauge in the SOV representation 

Let us recall that the SOV coefficients 3,(T]r) and d(?7r) are specified only by the average value 
relations (13.151) required by the cyclicity of the SOV representation. Then we can fix the gauge 
in the SOV representation fixing a couple a(A) and d(A) of Laurent polynomial solutions of the 
following average relations: 

p p 
AiA) + 7i3(A) = n ' ^(^) + '^^(^) = n ^(^'^) ' (3.22) 

k=l k=l 

where 7 and 5 are constants to be fixed. As it will be clear in the next sections, it is important 
to require that these solutions satisfy the further conditions: 

(a(A))* = d(A*), (3.23) 



13 

if aoE Za(A) ^ q'ao^ Za(A) VA;G {1,...,2/}, (3.24) 

Za(A)n Zj^i-i^^^^h^ = iD, (3.25) 

where Z/(a) denotes the set of the zeros of the function /(A). Let us denote with: 

Sp,Ao = (Ao, gAo, g^'Ao) E ap-string of center Aq, (3.26) 

then condition ( 13.241) . in particular, implies that a(A) and d(A) are free from p-strings0. Note 
that in the untwisted representations (f/„ = Ki = 1), we have defined these Laurent polynomials 
as solutions of (13.221) with 7 = — 5 = 1 which also satisfy (|3.23l) - (|3.25l) . This was possible as 
for untwisted representations it holds: 

^(A) = V{A), B{K) = C{A). (3.27) 

It is worth to note that this properties together with (13.111 ) imply that for untwisted representa- 
tions the average values of the elements of the monodromy matrix are real Laurent polynomial 
in A. However, this reality condition is lost for the twisted representations here considered; 
indeed, we have the following: 

Lemma 5. Almost for all the values of the parameters Kn, X„, Un, Vn the Laurent polynomials 
A{Iy) and T>{K) are not identical as well as -B(A) and C(A). 

Proof. The Lemma is true for N = 1, as it is clear from the one-site average formula (|3.21i 
Then, for N > 1, the proof follows by induction using Lemma |3] □ 

The previous Lemma and the Hermitian conjugation properties (13.1 II) imply that for generic 
twisted representations the only way to simultaneously satisfy (13.221 ) and (13.231 ) is by imposing 
7 = 5 = 0. Then representing: 

N 

^n(A) = \{{NZ^ - Z^jk) , (3.28) 

n=l 

we can chose to define d(A) by (13.231) and a(A) as any even Laurent polynomial of the form: 

N 

a(A) = aN n(^/< - with a^ = and (O^ = Z;f , (3.29) 

n=l 

with asymptotic s for N even: 

lim A±^a(A)= (3.30) 

logA->=Foo \ J-j^ / 

where 6* is a fixed p-root of the average value of the 6-charge. Moreover, we impose the 
following prescriptions: 



'°Here and in the following we say that a function /(A) is free from p-strings meaning that the set of its zeros 
does not contain p-strings, i.e. Sp_Ao ^/(a) ^-^o G C. 
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a) In the case of real or imaginary zeros of ^(A): 

lf3ne {1, ...,N} : Z;f G Z^(a), {Z^f eR ^ z^^e Z,(a), (z:)Vg G M. 

(3.31) 

b) In the case of multiple zeros of ^(A) with multiplicity r: 

If 3ni,...,n, G {1,...,N} : Z^^ = ... = Z^ G Z^(a) ^ = ... = G Z,(,). 

(3.32) 

c) In the case of couples of complex conjugate zeros of ^(A): 

If 3n^mG{l,...,N}:(Z;f)* = Z;^G Z^(A) -> {zl/ q)* = z^, e Z^^^y 

(3.33) 

Then under these conditions a(A) and d(A) satisfy the requirements (I3.231 )- (I3.251 ). 

3.3 Constructive proof of the existence of cyclic SOV representations 

In the following subsections we will show by recursive construction the 

Theorem 1. Almost for all the values of the parameters Kn, Xn,Un,Vn there exists a SOV 
representation for the SG model, in fact, we can prove that the one-parameter operator family 
B(A) is diagonalizable with simple spectrum. 

3.3.1 Recursive construction of B-eigenstates 

We will construct the eigenstates ( ?7 1 of B(A) = Bn(A) recursively by induction on N. 

In the case N = 1 we may simply take ( r/i | = {vi\, where ( f i | is an eigenstate of the operator 
vi with eigenvalue v. It is useful to note that the inhomogeneity parameter determines the subset 
of C on which the variable r]i lives, 7]i G {^i/vi)E>p. 

Let us assume we have constructed the BM-eigenbasis for any M < N, then the eigenstates ( r] \ 
of Bn(A) may be constructed in the following form 

^^1 = EE^n(^Ixz;Xx)(xJ®(xJ, (3.34) 

where ( Xz I and ( Xi I are eigenstates of Bm(A) and Bn_m(-^) with eigenvalues parameterized as 
in by the tuples Xi = {Xza)a=i,...,M and Xi = (Xia)a=i,...,N-M' respectively. It suffices to 
consider the cases where N — M is odd. 



\3.3.1\ a Dependence of the kernel K^{rj \XijX±) w.r.t. Xi <^nd Xi- 
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From the formula 

1 N — M 

(A) (3.35) 
it follows that the matrix elements of the kernel K^{ri \XxiXi_) have to satisfy the relations 

(A)B 1 N— M 

(A) + B,m(A)D 

1 N— M V I Xz, Xi ) 

N (3 36) 

= n T n (^/^'^ - ^"/^) ^ I X.; Xx ) , 

n=l a=l 

where we used the notation 0* for the transpose of an operator 0. Let us assume that 

XiaQ^' i ZdetqM,,N-M(A), XM^'' ^ ZdetqM,,M(A) and x±a(t^ ^ Xihq^\ (3-37) 



where hi G a G {1,...,N — M} and h G {1,...,M}. Under these assumptions the 

previous equations yield recursion relations for the dependence of the kernel in the variables 
Xia and Xzb simply by setting A = Xia and A = Xzfe- Indeed for A = Xia the first term on the left 
of (|3.36l) vanishes leading to 



N-M . M] 

'^laKAv\Xz\X^) dM~^xJ Xm^ n —W^xJXzb-Xzb/x^ 



1 '^ri 1 
n=l a=l 



/n 



(3.38) 



{xJ%-%/x 



6=1 



= K^{v\Xz;xJ vl 

while for A = Xia one finds similarly 



M . N-M 



1 ^'i !, 1 

n=l b=l 



(3.39) 



iXza/Vb-Vb/X 



b=l 



where d-^(Xia) and a^{Xza) are the known coefficients of the SOV representations in the sub- 
chains 1 and 2. If M is even we find the recursion relation determining the dependence on Xzm 
by sending A — )■ oo in (13.361) . leading to 



M-l , N-M 



tL^.(^Ix.;xJ— n— n — = K,iv\Xz;v.)ll-- (3.40) 

XzA a=l 6=1 -^1^ 6=1 



\3.3.1\ b Dependence of the kernel K^{rj \XijX±) w.r.t. rj. 
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The action 

{r]\Dt,{r]i) = -detqM^ M(g?7i)(?7|BjJ^(g?7i)Bi n-mI??*)) 
follows from the formula 

1 N — M 1 N— M (A) 



(3.41) 



(3.42) 



when we express m(A) by the quantum determinant detqMi m(A) in the subchain i and we 
use that 

{r]\B^ j^-uiVi) = -(^|Di j,-M{Vi)KMiQVi)^i m{qVi)- (3.43) 

So the kernel matrix elements K^{ri | Xz! Xi ) have to satisfy the following recursion relations 
for the dependence w.r.t. rj 



V I X2; Xi )detqM^ ^{qr]i) Yf^J^it^ali) {rji/Xia - X^a/v 



XIm Ua=N-M+l l^al i ]\^^^^ {qVi/Xia " Xia/qV: 



[M] 



(3.44) 

If N is even we have to consider further the asymptotic of the equation (13.421 ) which by the 
asymptotic behavior, analyzed in appendix |Bj leads to the identity: 

lim A±^D,(A) = - lim A^^detqM, M(Ag)B-L(Ag)B, ,_m(A). (3.45) 

log A— ^=poo log A— >=poo 

From which we get the following recursion relation determining the kernel dependence on r^^: 



N-l 



M ^2 N-M 

Tn^n( V I Xz; Xi ) n = ^n( v\Xz;Vi) Yl JJ ^ 



a=l 



b=l 



X.b 



1 b ■ 



(3.46) 



b=l 



3.3.2 Determination of gauge-invariant SOV dates: Zr, Z/\, Zd, A^{Zr), T>^{Zr 



The condition of compatibility among the recursion relations (I3.38I) - (I3.39I ) and (13.401) and the 
requirement of cyclicity, (T^^)^ = (T^^)^ = 1, provide a system of N algebraic equations in the 
N unknown Za = Va' 



N-M 



M] 



^^xn-m(xL) ixlZr n ^U^^'Jx'lb-xlJx'J = ivTrU^^Uv'b-vUx'J, (3.47) 



1 I. 1 
n=l 6=1 



b=l 



M .„ N-M 



Auixla) n - n (^'a/^'b - XlbMa) = ivTY U^^UVI - Vllxla) , (3-48) 

1 '^'^ L 1 11 

n=l b=l b=l 

with for M even 

^ M-l -| N-M ^ N ^ 

" " " (3.49) 



XiA a=l 6=1 % 
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The simplicity of the spectrum of B(A) in the SOV representations of both the subchains i and 
2 implies that the above system of equations completely determines the unknown Za in terms of 
Xza' X?a- Indeed, we can reformulate this system of equations in terms of the following Laurent 
polynomial equation: 

(A/Z,-Z,/A), (3.50) 

n=l " a=l 

where the l.h.s. is formed out of the known average values of the monodromy matrix elements in 
the subchains i and 2. Therefore the problem to determine the unknown for a G {1, [N]} 
is reduced to the problem to determine the zeros of the known Laurent polynomial at the l.h.s. 
of (13301). 

The requirement of cyclicity = 1 determines the remaining gauge-invariant SOV dates: 

V^{Z,) = -det^Mum^^^^, ^NiZ.) = -det^M^.uiZ,) ^""^f^ (3.51) 

for any i G {1, [N]} where we have defined: 

p 

det^Mx{A) = Y[det^Mx{q''X), X = M,N-M,N, (3.52) 

a=l 

while for N even the condition (Ti^)^ = 1 reads: 

N-l N-M N -^2 

a=l 6=1 b=l 

Note that from the identity: 

detq7Wx(A) = Ax{^)Vx{^) - Sx(A)Cx(A), X = M, N - M, (3.54) 

the equations (13.511) and (|3.53l) are equivalent to the equations (13.191 ) and (13.201) as well as the 
equation (13.501) coincides with (|3.17l) of Lemma [H 

3.3.3 Completeness of B-eigenstates and simplicity of B-spectrum 

Here we show that the set of B-eigenstates ( rj \ constructed in the previous subsection is com- 
plete, i.e. it defines a basis of the representation; this is done showing that there are distinct 
corresponding B-eigenvalues. In particular, we can prove: 

Proposition 2. The SOV dates Zr with r G {1, [N]} are all distinct for almost all the values 
of the parameters Kn, X„, f/„, Vn of the SG model. 



^M(A)i3N-M(A) + Sm(A)I)n-m(A) = Z^- TT 
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Proof. Note that we can apply Lemma |4] to the average values on the l.h.s. of (13.501) and so 
we have that the r.h.s. of (|3.50l) times Y^=i ^ Laurent polynomial of maximal degree 

1 in each of the parameters Kn, X„, f/„, V„. The consequent functional dependence^ of the 
Z^, . . . , Zjj^j w.r.t. these parameters implies that it is sufficient to show the non- vanishing of the 
Jacobian: 

J{X- K, U, V) = det ^ (3.55) 

\ s/ r,s=l,...,[N] 

for some special values of the parameters K^, Un, Kt in order to prove that J{X; K,U,V) ^ 
for almost all the values. Whenever J{X; K, U, V) ^ 0, we have invertibility of the map 
Z = Z(Xi, . . . , X[N]) from which the claim of the Proposition follows. 

To show that (|3.55l) is indeed satisfied, let us choose Kn/Vn = for n = 1, [N], then the 
average values (13.211) of the Lax operators simplify to 

C(A) = ±( ° V^iAVjX^-XjV^A) 

V K(A/KX„-X„K/A) 

Using now (|3.12l) to compute the l.h.s. of (|3.50l) we get: 

[N] 



n=l a=l 



{A/Za-ZJA), 



(3.57) 

from which the fact that J(X; K, U,V)^{) trivially follows. □ 

We conclude that our construction of B(A)-eigenstates will work if the representations TZ^, TZm 
and 7?.N-M are all non-degenerate. Theorem 1 follows by induction. 



4. Characterization of T-spectrum 
4.1 SOV characterization of T-spectrum 

In the SOV representations the spectral problem for T(A) is reduced to the following discrete 
system of Baxter-like equations in the wave-function \l't(?7) = ( 1 1 ) of a T-eigenstate \ t): 

t(r7,)v[/(r/) = a(r/,)vI/(T;(r/)) + d(r7,)v[/(T+(r7)) Vr G {1, [N]}, (4.1) 

where (r/i, //n) G Bn with Bn the set of zeros of the B-operator in the SOV representation. 
Here we have denoted with Tf{rf) = (r/i, . . . , q^^rjj., . . . , //n) and a(?7r) and d{rir) the coeffi- 
cients of the SOV representation as defined in subsection l3.2.2[ In the case of N even, we have 

11 Let CT.L^1(Z) be the de gree n elementary symmetric polynomial in the variables Z^, then ctL^' (Z)/ CTjJJj (Z) are 
Laurent polynomials of maximal degree 1 in each one of the parameters Kn, Xn, Un,Vn- 



to add to the system (|4.1I) the following equation in the variable ?7n: 
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T^^kiv) = q-'^^kiv), (4.2) 

for t{\) e Sy^ with k E {0, 21} and 9 = rjA Y^Zl Va/ Y^=i ^a, obtained from the asymp- 
totics of T(A) given in appendix |Bl 



4.2 Characterization of T-eigenvalues as solutions of a functional equation 



Let us introduce the one parameter family D{X) of p x p matrix: 

/ t{X) -d(A) ••• 

-a(gA) t{qX) -d(gA) 







\-d(g2'A) 



-a(A) \ 






-a(g2'-iA) t(g2'-iA) -d{q^'-'X) 
-a(g2^A) tiq^'X) J 



(4.3) 



where for now t{X) is just a real even Laurent polynomial of degree N in A. 

Lemma 6. The determinant of the matrix -D(A) is an even Laurent polynomial of maximal 
degree N in A = A^. 



Proof. Let us start observing that D{Xq) is obtained by -D(A) exchanging the first and p-th 
column and after the first and p-th row, so that 

det D{Xq) = detD{X) VA G C, (4.4) 
p p 

Being a(A), d(A) and t{X) even Laurent polynomial of degree N in A then the parity of 
detpD{A) and the fact that it is a Laurent polynomial of maximal degree N trivially follow. □ 



The interest toward the function detp -D(A) is due to the following: 
Lemma 7. Let t{X) G Sj, then t{X) is a solution of the functional equation: 



detD(A) = 0. 
p 



(4.5) 
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Proof. Note that from t{\) G Sj and the SOV characterization of the T-spectrum we have that 
the requirement that the system of equations (|4.1I) admits a non-zero solution reads: 

deiD{ril) = Q VaG {1,...,[N]} anrf (r/i,...,r/[N]) GBn, (4.6) 
p 

In the case of N even, we have moreover: 

Um A^^detL)(A) = 0, (4.7) 

logA— >=Fco p 

which simply follows by observing that: 

lim A±^detD(A) = det + g-^^^^^^^^^.j+i - {qH + q-''e-^))5iA\ 

log A— >=poo p p " ' ' 



N 

X ^ ^] 

n=l 



1) n = 0- (4.8) 



for t{X) G St''' and A; G {0, 2/}. 

Then the function detpD(A) is zero in N+cn different values of A^ which thanks to Lemma [6] 
implies the statement of the Lemma. □ 

Remark 2. Let us note that the same kind of functional equation det D(A) = also appears 
in [|48ll49l[50l[5n . There it recasts, for different integrable quantum models at the roots of unit, 
the functional relations which result from the truncated fusions of transfer matrix eigenvalues; 
this is in particular true for the r2-mode0 ESllSTllSll 153)1 . 

4.3 Construction of Baxter equation solutions from T-eigenvalues 

Theorem 2. For any t{\) G Sj, we can construct uniquely up to normalization a real polyno- 
mial: 

2lN-bt 

QtW = n - ^')' < < 2/, < 6i < 2/N, (4.9) 

h=l 

which is a solution of the Baxter functional equation: 

t{X)Qt{X) = a(A)Qi(Ag-i) + d{X)Qt{Xq) VA G C, (4.10) 
free from p-strings. Moreover, for N even and t(A) G Sj'^ with k G {0, 21}, it holds: 



at = k, bt = k mod p. 
The SOV representations of this model were analyzed in a series of works ll54l l55l l56l . 



(4.11) 
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Proof. Construction: Let us notice that the condition t{\) E Sj implies that the p x p matrix 
D{X) has rank 21 for any A G C\{0}. Let us denote with 

C„(A) = (-l)^+MetA,,(A) (4.12) 

the cof actor of the matrix -D(A); then the matrix formed out of these cofactors has rank 1, 
i.e. all the vectors: 

V,(A) = (Q,i(A), Q,2(A), C,,2i+i{\)y G V2 G {1, 2/ + 1} (4.13) 
are proportional: 

V,(A)/Q,i(A) = V,(A)/C,,i(A) V2,j G {1, 2/ + 1}, VA G C. (4.14) 

The proportionality (14.141 ) of the eigenvectors Vj(A) implies: 

C2,2(A)/C2,i(A) = Ci,2(A)/Ci,i(A) (4.15) 

which, by using the property (|A.2I ). can be rewritten as: 

Ci,i(Ag)/Ci,2i+i(Ag) = Ci,2(A)/Ci,i(A). (4.16) 

Moreover, the first element in the vectorial condition D{\)Ni{\) =0 reads: 

t(A)Ci,i(A) = a(A)Ci,2z+i(A) + d(A)Ci,2(A). (4.17) 

In appendix lAl we have analyzed the properties of these cofactors and thanks to Lemma |9] and 
[10]we can introduce now the polynomials Ci^iCi_i(A),Ci 2Ci,2«+i(A) and Ci,2Ci,2(A) defined by 
simplifying the common factors in Ci i(A), Ci 2/+i(A) and Ci 2(A), respectively. In particular, 
from Lemma |9] and [10] we have that Ci i(A) is an even polynomial of the form: 

Ni,i 

Ci,i(A) = n(A' - AD, Ni,i < 2/N, (4.18) 

h=l 

which furthermore satisfies the properties: 

Reality : (Ci,i(A))* = Ci,i(A*) and i Zc^^(;,), S^.Ao 'Z^c,,,(x)^ VAq G C. (4.19) 
Now, being by definition J ^(;,)nZc^^( A) = Zc^^(A)nZc^ 2,+i(a) = 0, equation dUJ) implies: 

Ci,2m(A) = g'^^'^Ci,i(Ag-i), Ci,2(A) = g-2^"Ci,i(Ag) , ip = C^,^/c^,2 aphase, 

(4.20) 

and then equation (|4.17l) assumes the form of a Baxter equation in the polynomial Ci_i(A): 

t(A)Ci,i(A) = a(A)Ci,i(Ag-i) + d(A)Ci,i(Ag), (4.21) 
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with coefficients a(A) = (g2Ni.iy,)a(A) and d(A) = (g2Ni,i^)-id(A). 

Let us denote with -D(A) the matrix defined as in (14.31) but with coefficients a(A) and d(A). 
Then the consistence condition detp Z)(A) = of the Baxter equation (14.211 ) and the condition 
t(A) e Ej implies: 

det D{A) - det ^(A) = < — > {^p" - 1) U(A) - ^^"^©(A)) = 0, (4.22) 
p p 

which is true if and only if is a p-root of the unity. Note that (14.221) is derived by using the 
following expansion for detp -D(A) and a similar one for detp Z)(A): 

det D{A) = AiA) + ViA) - a(A)d(A/g) det /^(i,2i+i),(i,2«+i) (A) 

-a(Ag)d(A) det Z}(i,2),(i,2)(A) +t(A)det A,i(A), (4.23) 

where D(^h,k),{h,k){^) denotes the (2/ — 1) x (21 — 1) sub-matrix of -D(A) obtained removing the 
rows and columns h and k, plus the formulae^: 

det 1 ( A) = det Di 1 ( A) , (4.24) 

det Ai,2),{i,2)(A) = det D(i,2),(i,2)(A), (4.25) 
det I?(i,2i+i),(i,2«+i) (A) = det £'(i,2«+i),(i,2«+i) (A) . (4.26) 

Finally, we can define our polynomial solution (14.91 ) of the Baxter equation (14.101 ) by: 

Qi(A) = A'^'Ci,i(A), (4.27) 
where g""* = g^Ni,!^^ ^ {g^ ..^ 21} and h = 2/N - Ni,i. 

In the case N even, for t(A) G with k G {0, 21}, we get from the Baxter equation and 
the asymptotics of the coefficients (13.301) the following asymptotic conditions: 

„,„«£[^ = ,^ „,„ «iM = (4.28) 

Then the characterization (14.111 ) is just a corollary of (14.281 ). □ 

Proof. Uniqueness: Let t(A) G Sj and let Qt{\) be the solution of the Baxter equation (14.101) 
constructed in the previous part of the proof. Then denoting with Qt{X) G C[A] any other 
solution, we can define the q-Wronskian: 

Wt{X) = Qt{X)Qt{q-'X) - Qt{X)Qt{q-'X) ■ (4.29) 



'^They follow from Lemma 3 of |2 | thanks to the tridiagonality of these matrices. 



which by the Baxter equation satisfies the equation: 
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(4.30) 



Thanks to the cyclicity (T+)p = 1 the average of the above equation reads: 
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{A{A) - V{A)) WtiA) = with Wt{A) = J] WtiXq''), 



(4.31) 



which by Lemma[5]implies W{\) = 0. It is then easy to see that this implies that (A) = Qt(A) 



Remark 3. The previous theorem implies that for any t(A) G Sj the polynomial solution 
(5t(A) of the Baxter equation can be related to the determinant of a tridiagonal matrix of finite 
size p — 1 for the rational = p' /p. It is then relevant to remark that others determinant rep- 
resentations for Baxter equation solutions are known in literature; for example in the quantum 
periodic Toda chain. There, these solutions are expressed in terms of linear combinations of 
determinants of semi-infinite tridiagonal matrices [|26l l57l[58l . The above remark explains the 
interest in a careful analysis for the SG model of the limit — )■ with irrational. Indeed, 
under this limit (i.e. p', p — )■ +oo) the dimension of the representation in any quantum site 
diverges as well as the size of the tridiagonal matrix associated to a solution (5t(A) of the Baxter 
equation. Then, it is natural to investigate if characterizations like those encountered for the 
quantum periodic Toda chain apply to the SG model for irrational too. In particular, this is 
interesting for a potential reformulation0 of the SG spectrum in terms of solutions of nonlinear 
integral equations (NLIE) as shown recently in [i6ri for the quantum Toda chain. On the other 
hand, it is worth remarking that the SG spectrum admits another NLIE reformulation which 
IS of DDV-type0 and satisfies the feature of completeness as consequence of the Bethe ansatz 
completeness proven in Subsection 14.61 

4.4 Construction of T-eigenstates from T-eigenvalues 

The results of the previous subsections allow both to introduce a complete characterization of 
the set St and to construct one T-eigenstate \t) for any t{\) E Sy; indeed, we can prove the 
following: 

'"^The use of NLIE to reformulate the spectrum characterization in integrable quantum models goes back to 
Il59l l60 l. NLIE reformulations similar to that presented in [61 1 also appear in ll TSl l62l l45l l63l . 

This type of NLIE are derived by a reformulation of the Bethe ansatz equations and they were introduced and 
analyzed for fermionic lattice regularizations of the Sine-Gordon model in Il64l - ||7T1 and in ll72l - ll75l for a related 
model. 



up to normalization when we have assumed Qt{^) free of p-strings. 



□ 
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Theorem 3. The set Sj coincides with the set of all the t{X) G M[A^, A^^]n/2 solutions of the 
functional equation: 

detD(A) = 0, VAgC. (4.32) 

p 

Moreover, for any t{X) E T,j a corresponding T-eigenstate is characterized by: 

N 

^t{r]) = {r]i,...,7]^\t) = Y[Qt{r]r), if N zs odd, (4.33) 

r=l 

N-1 

^t{v) = {Vi,-,VN\t)=v^''Y[Qt{Vr), tf N ts even and t{X) e (4.34) 

r=l 

w/iere QtiX) is the Baxter equation solution constructed in Theorem^ 

Proof. The inclusion of Sj in that set of solutions of (14.321 ) was proven in Lemma |7J To prove 
the reverse inclusion of sets, we have just to observe that the construction of Qt(A), presented 
in Theorem m holds for any t{X) G M[A^, A"^]n/2 solution of (|4.32l) . Then, the t(A) and the 
\E'( (77), (14.331) and (14.341) . are solutions of the system of Baxter like equations (14.11) and (14.21) and 
so they define a T-eigenvalue and a corresponding T-eigenstate. □ 

4.5 Simplicity of T-spectrum 

In this section we show that the spectrum of the transfer matrix T(A) is non-degenerate (or 
simple). Let us start proving the following: 

Lemma 8. For almost all the values of the parameters Kn, Xn, Un, Vn, the average values of 
the monodromy matrix elements A(A) and D(A) satisfy the inequalities: 

A{Za) ^ V{Za), Va G {1, [N]}, (4.35) 

where are the zeros of the average value o/B(A). 

Proof. We have to prove that the functions: 

7a[K^, X„, t/„, K) = ^N(^a) - V^[Za), Va G {1, [N]} (4.36) 

are nonzero for almost all the values of the parameters Kn, X„, Un,Vn. The Lemma |4] and the 
functional dependence^ of the Z^, . . . , Z^-^-^ w.r.t. these parameters implies that it is sufficient 
to show that the functions J-'a are nonzero for some special value of the parameters in order to 
prove our Lemma. Note that the following identities hold: 

J^a \u„ = l Vne{2,...,N},y„ = l Vne{l,...,N}= (1 + ~ f^f^ )^2, N-1 (^a) , (4.37) 

"■See footnote [TTI 
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Here, we have used the decomposition of the chain in a first subchain i, formed by the site 1, 
and a second subchain 2, formed by the remaining sites. Moreover, we have used the identities: 

(A,i-^^i,i) = {l + Kl){U^-U^'), B,,i(A) = C,,i(A), (4.38) 

S,,N-i(A) = C,,N-i(A), P,,N-i(A) = A,N-i(A), (4.39) 

which follows by using Lemma [3] and (13.211) for the special choice of the parameters done in 
(14.371) . Then for Ki ^ ±i and U\ 7^ ±1, we have only to show that we can always take 
A,N-i(^a) ^ for any a e {1, [N]}. 

Note that from the identities: 

detq7W,,N-i(A) = (A,N-l(A))'-(i3,,N-l(A))^ (4.40) 
i3N(A) = A,i^.,N-i(A) + S,,i(A)A,N-i(A), (4.41) 

we have that ^2,N-i(^a) = if and only if Za is a double zero of detqA^^ n-i(A). However, 
this is not the case in general because the function detqA^2,N-i(A) has not double zeros for 
general values of the parameters X„ and K^, as it is clear from the formula: 

N 

detqA^,,N-i(A) = n n ^ni^/^n,h " M„,;,/A), Mn,± = ±^K^^Xr,. (4.42) 

h=±l n=2 

which is obtained averaging the quantum determinant (|3.61) . □ 

Theorem 4. For almost all the values of the parameters Kn, Xn, Un, Vn of a twisted represen- 
tation, the spectrum o/T(A) is simple. 

Proof. We have to prove that up to normalization for any given t(A) G Sj the wave-function 
(14.331) and (14.341) are the only solutions of the system (143]) and (|43l) . 

Let us denote with ^tij])^ 77 G Bn, any other solution corresponding to the same T-eigenvalue 
t(A). Then, we can define the q-Wronskian: 

WtM = ^t{T;{r])) - ^t{r]) ^t{T;{r])), Vr G {1, [N]}. (4.43) 

which by the system of Baxter like equations (14.11) satisfy the equations: 

a(r7,) WtAv) = ^iVr) WtA'^tiv)), Vr G {1, [N]}. (4.44) 

Thanks to the cyclicity (T+)p = 1, the averages of the above equations read: 

21 

{A{Zr) - ViZr)) WtAv) = with WtAv) = n ^tA'^riv)), (4.45) 

A:=0 
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which by Lemma [8] implies: 

W,AV) = ^ ^^J^ = ^ ®N and Vr G {1, [N]}. (4.46) 

In the case N odd, these identities just imply the uniqueness of the solution of the system (|4.1I) 
and its factorized form given in (I4.33I ). In the case N even, they imply the factorization: 

N-l 

^t{v) = ft{VN)l[Qt{Vr). (4.47) 

r=l 

Then we get /t (?7n) oc t]^^ as the unique solution of the equation (14.21) in t^n for t(A) G Sj'^. □ 

Remark 4. It is worth to point out that for the twisted representations of the SG model the 
spectrum of the transfer matrix is simple both for N odd and even. This is due to the fact 
that by definition for N even these representations are characterized by an average value of the 
6-charge different from 1, i.e. 9^ ^ 1. 



4.6 Completeness of the Bethe ansatz 

An important consequence of the previous analysis is that it naturally leads to the complete 
characterization of the transfer matrix spectrum (eigenvalue and eigensate) in terms of real 
polynomial solutions of the associated Baxter functional equation. 



ies that to any t(A) G Sj it is associated one and only 



Let us observe that Theorem |2] imp 

one self-complex-conjugate solutior[|J (±Ai, . . . , ±A2iN_fe) of the following system of Bethe 
equations: 

Conversely, let us consider a self-complex-conjugate solution of the Bethe system of equations 
(14.481) . then it defines uniquely a real polynomial Q{\) by the equation (14.91) . Now, by using 
Q{\) we can construct the function: 

t(A) = (a(A)Q(g-iA) + d(A)g(gA))/g(A) . (4.49) 

which, thanks to the Bethe equations (14.481) . is nonsingular for A = ±\k, Vfc G {1, . . . , 2/N— 6} 
and is a real even Laurent polynomial of degree N in A. Moreover, we can also uniquely 
construct a state 1 1 ) by inserting Q{\) in equation (14.331 ) for N odd (and (14.341 ) for N even). 
Then, t{\) and the wave-function '^tiv) = ivi^ ■■■^Vn It) satisfies by definition the system of 
Baxter like equations (14.11 ) (and (14.21 ) for N even). This implies that starting from the given 



'^It means that the tuple (±Ai , . . . , ±A2/N-b) is invariant up to reordering under complex conjugation. 
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solution of the Bethe equation we have uniquely reconstructed the T-eigenvalue t{\) and the 
corresponding T-eigenstate \t). This establishes a one-to-one correspondence between the 
above solutions to (|4.48l) and the spectrum of the transfer matrix (Completeness of the Bethe 
ansatz). 



4.7 Reconstruction of Q-operator 

In the previous sections, we have given a complete and operative characterization of the spec- 
trum of the transfer matrix T(A) starting from its SOV characterization. It is worth pointing 
out that we have derived our results without any need to introduce a Q-operator. Anyhow, it is 
worth noticing that our construction has as built-in property the existence and characterization 
of the Q-operator. This feature is even more relevant in the case of twisted representations of 
the SG model where a direct construction by cyclic dilogarithm functions can be applied only 
on the sub-variety of representations generated by points on algebraic curves. 



Definition 1. Let Q(A) be the operator family defined by 

Q{X)\t) = Qt{X)\t), where Qt{X) = X^^'-'^-^^QtiX), Vt(A) G Ey. (4.50) 

Here, Qt(A) is the real polynomial corresponding to t(A) by the injection defined in Theorem^ 
and \t) is the corresponding T-eigenstate. 

Then the following theorem holds: 

Theorem 5. The operator family Q(A) is a Baxter Q-operator: 

(A) Q(A) and T(A) satisfy the commutation relations: 

[Q(A), T(^)] = [Q(A), Q(/i)] = VA, ^ e C, (4.51) 

plus the Baxter equation: 

T(A)Q(A) = a(A)Q(Ag-i) + d(A)Q(Ag) VA e C, (4.52) 

and for N even: 

[Q(A),e] = VAgC. (4.53) 

(B) Q(A) is a polynomial of degree 2/(N + 1) in X^: 

2«(N+1) 

Q(A)= J2 Q"^'"' (4.54) 

n=0 

with coefficients Q„ self-adjoint operators. 



'^Here, we have decided to define Qt{X) by multiplying the function Qf (A) for the p-string when at is odd. 
This is done to get that Qt(A) is always a real even polynomial in A which is still a solution of ( 14.101 ). 
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Proof. Note that the self-adjointness of the transfer matrix T(A) implies that Q(A) is well de- 
fined being its action defined on a basis. The commutation relations and the Baxter equation 
are trivial consequences of Definition 1, thanks to the simplicity of the T-spectrum. Finally, the 
reality condition of the Q-eigenvalues: 

{Qt{\)y = Qt{\*) VAgC (4.55) 
implies that the operators Q„ are self-adjoint. □ 

Remark 5. In ll53ll . Baxter has extended the construction of the Q-operator by gauge trans- 
formations of [|25l to the r2-model for general cyclic representations not restricted to those 
parameterized by points on algebraic curves. The main tool there was a generalization of the 
discrete dilogarithm functions. In particular, he has remarked that asking the cyclicity only for 
the products of couples of these functions the Q-operator can be still well defined for general 
cyclic representations of the model. Note that there are no doubts that this construction can 
be adapted to our general representations of the SG model leading to the representation of the 
Q-operator in the basis (|2.5I ) of the representation. However, it is also clear that this last charac- 
terization of Q cannot add any information to our construction of the transfer matrix spectrum 
but it can be potentially useful for more physical aims. In particular, we can try to use it to 
investigate the lattice dynamics [|76l iTTll associated to the general cyclic representations of the 
SG model. This can be done following the the same line drawn in [IJ, i.e. showing that the 
discrete evolution operator can be written in terms of the Q-operator. Here, the main physical 
aim is to prove that these representations of the SG model define lattice regularizations of the 
Sine-Gordon model for general a-sectors of the quantum theory (see Section 2.2 of fT] for the 
definition). 

A. Properties of the cof actors Cj j(A) 

Let t{\) be just a real even Laurent polynomial of degree N in A, then the cofactors of the matrix 
D{X) satisfy the following: 

Lemma 9. The cofactors Cjj (A) are even Laurent polynomials of maximal degre^ 21N: 

h=l 

which satisfy the following properties: 

Ch+i,k+i{\) = Ch,k{W) V^, h,ke {1, (A.2) 
'^Theaij- and 6^ j- are nonnegative integers and by definition A(ij)ji ^ Oforany/i G {1, 2ZN— (0^^ +6^^)}. 
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and 

(Ci,i(A))* ^ Ci,i(A*), (Ci,2(A))* ^ Ci,2,+i(A*). (A.3) 

Proof. The characterization (lA.ll ) of the cofactors trivially follows from being a(A), d(A) and 
t(A) even Laurent polynomial of degree N in A. Note that by the definition (14.121 ) of the cofac- 
tors Cjj(A) the equations (IA.2I) are simple consequences of = 1 and are proven exchanging 
rows and columns in the determinants. 

Now let us prove the property (|A.3I ) for the cofactor Ci,i(A) = det2/-Di,i(A), where: 

Z}i,i(A) = \\t{\q^')6H,k - 3.{Xq^)6h,k+i - ^iM'')Kk-i\\^<^h<2i,i<k<2i ' 
then by the properties under complex conjugation of t(A), a(A) and d(A) it holds: 

(Z}i,i(A))* = ||t(A*g'^)5,,, - d(A*g^-^)5,,fe+i - ^{>^*q'-'')Kk-i\\,^j,^^i^,^^^,i , (A.5) 
Let us denote with the 21 x 21 matrix of columns: 

VaG{l,...,2/}, (A.6) 

where stays for the column a of the matrix X, and similarly with D^'^ the 21 x 2/ matrix 
of rows: 

Ra''^ = R^l'^, Va e {1, 21}, (A.7) 
where i?^ stays for the row a of the matrix X. Then we get the identity: 

D^f (A*) ^ (Di,i(A))* ^ (Ci,i(A))*^Ci,i(A*), (A.8) 

being the determinant invariant under an even number of column and row exchanges. In a 
completely similar way, we can prove the identity: 

(det D(i,2).(i,2)(A))* = det Z^(i,2),(i,2)(A7g), (A.9) 
then by using the cofactor expansions: 

21 



Ci,2m(A) = n^(^^') + ^(Vg)det^(i,2),(i,2)(A/g), (A.IO) 

h=l 
21 

Ci,2(A) = J]d(Ag") + a(Ag)detD(i,2),(i,2)(A). (A.ll) 

h=l 

the second identity in (IA.3I ) is a corollary of (IA.91 ). □ 
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Lemma 10. Let t{X) G Sj, then the following identities hold: 

0-1,1 = ai,2, &i,i = bi,2 (A.12) 
where aij and bij are the nonnegative integers defined in i4.iD . 

2ci,i(A) n Zci_2(A) = Zci,i(A) n Zc,_2,_^,(A), (A.13) 

3 Sp,x, C Zc,,,(A) ^ ^p,Ao nZc,,2(A) ^ 0, (A.14) 
/or ^p^Ao = (Ao, gAo, g^'Ao) any p-string. 

Proof. Under the assumption t (A) G Sj we have proven the validity of the equations (14.161 ) and 
(14.171) . Then the identity (IA.12I) follows by using (14.161) and the property (IA.3I ). Concerning the 
proof of identity (IA.13I ) thanks to the prescriptions (|3.231 )- (I3.25I) . chosen in subsection |3.2.2l we 
can follow word by word the proof given in Lemma 5 of the paper [2J. 

Let us assume by absurd that (|A.14I ) is not satisfied, then (14.161) implies that S^^Aq C Zc^ 2i+i(A)- 
Now this last condition implies that (14.171 ) can be satisfied only if Sp Ao C Z(i(A), which does not 
hold by the definition of d(A). □ 



B. Asymptotics of Yang-Baxter generators 

From the known form of the Lax operator we derive the following asymptotics for A — )■ +oo 
and of the generators of the Yang-Baxter algebras. 

N odd: The leading operators are Bn(A) and Cn(A) with asymptotics: 

(N \ / N (-1)1+" N \ 

n - n - ^'"^ n U-'^^ + ^ub-leadlng, (B. 1) 

a=l V V a=l a=l ) 

/ N \ / N (-1)'' N \ 

Cn(A) = n - n ^ - n + s^b-Mng. (B.2) 

\a=l V V a=l ^'^ a=l J 

N even: The leading operators are An(A) and Dn(A) with asymptotics: 

/ N \ / N (-1)" N \ 

An(A) = ( n y ) ( A^ n ^ + ^"'^ n ) + ^ub-Ming, (B.3) 



\a=l / \ a=l a=l 

(N \ / N (-1)1+" N 

n T n + ^'"^ n ^^^^i-'^^ ) + ^^b-leadlng. (B.4) 

a=l ) V a=l a=l 
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Note that these asymptotics imply for the SOV representation of the Yang-Baxter generators the 
following formulae§°: 



Nodd: 



N even: 



va=l 



N-1 



a=l 

N 



sov 



sov\-i o-i^sov 



, a=l J 
' N-1 \ 

n ^ 



a=l 



a=l 



Note that taking the average value of the last two formulae we get for N odd: 



a=l 



a=l 



while for N even: 



where 



N-1 N 



za = (6) n ^-^ n^- = ^A(e)- 



a=l a=l 
N 



1) 



1 + a 



a=l 



is the average value of the charge 0. 



(B.5) 



(B.6) 
(B.7) 



(B.8) 



(B.9) 



(B.IO) 



C. Comparison between the r2 -model and the SG model 

In this appendix we present the comparison between the SG model and the so-called r2 -model 
introduced and analyzed in the series of papers [|25l ISTI- llS&l . 



C.l Lax operators 



The Lax operator which describes the r2 -model reads: 



TT2 — A I = (c^n.+Xn + (i„^_X„^) Hn = [hn-'^n + ^n.+X^"*^) Y,„ i i 

Gn = (5'n,-X„ + (7r!,,+X„ j Y^ -Fn = {fn-^n + fn^+^n ) 



(C.l) 



^°Here, we denote with w^°^ the invertible matrix which defines the change of basis in the representation from 
the basis obtained by the tensor product of the v„-eigenbasis to the B-eigenbasis. 
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where: 



A. - r; ° 1 (C.2) 



and: 

9n,+ = fn,+dn-/hn,+ , Qn,- = fn,-dn,+/hn~- (C.3) 

Here, we have denoted with Y„ = Un/un and X„ = v„/f„ the generators of the local Weyl 
algebra Y^X^ = gX„Y„ with the following central values X^ =Y^ = 1. 

Then, the Lax operator of the r2 -model can be related to one of the SG model by: 

LI' = L'fa,, (C.4) 
when the parameters in (|C.ll ) are identified by: 

(C.5) 
(C.6) 
(C.7) 
(C.8) 

It is worth pointing out that in our SG model, we are asking that the elements of the Lax operator 
generate the algebra Uq{sl2) and that such a requirement implies the relations (IC.5I ). 

It is worth pointing out the different role played by the parameters m„ in these two models. From 
(IC.ll) . we see that the parameter m„ enters in the r2-Lax operator as a simple gauge while this is 
not the case for the SG Lax operator beinj^ [(Ti,A„] ^ 0. In particular, this means that while 
we are free to change uniformly along the chain the average values of the local generators u„ 
without any modification in the r2-transfer matrix this is not the case for the SG transfer matrix. 



d-n,— 




fn, 




fn,+ 




dn, 


+ = ~iXnVn, 


hn- 


= q ^ KnVm 




- q i^n '^n 


9n,+ 


— q f^nUn J 


9n- 


— q K„ Un- 



C.2 Map from the SG model to the r2-model for even chain 

The relation (|C.4I) directly implies that for an odd chains the transfer matrices in the two models 
have different spectrum. In the case of even chain a more precise analysis is need and it is here 
developed. 

Let us fix N = 2M with M positive integer; then, it is simple to define a map from the generators 
of the Yang-Baxter algebras in the SG model to those in the r2-model. We have just to point out 
that the relation (|C.1I ) implies: 

M"2(A) = (Tin,2 (M^^(A))(Ti ^ T"2(A) = n,, (T^^(A)) . (C.9) 

^'The presence of the cti tells us that what is a simple gauge for the T2-model is not a gauge for the SG model 
and vice versa. 
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where we have defined: 

2m— 1 \ '^2m—l ) 5 

t 

(C.IO) 

for m E {1, M} and a E {0, 1}; note that in the last equality we have used (I2.14I ). Then, we 
can characterize the map by its action on the generators of the local Weyl algebras: 

M 

= n^2n, where: FniUn) = u~\ FnK) = (C.U) 

n=l 

Until now, the discussion presented is common for both the untwisted and twisted representa- 
tions of the SG model. However, the distinction between these two cases appears evident when 
we observe that: 



Untwisted representations: The map is realized by the unitary transformation tt 



r2- 



M 



n,,(M^^(A)) = 7r,2M^^(A)7r,2, 7r,2 = JJfisn, (C.12) 

71=1 

after the flipping ,^2n-a — > (— e )^^~"^^2n-a of the inhomogeneities. Here, e is the sign 
defined in (12.121) and are the local unitary transformations defined by the following 
action: 

on the elements of the basis (12.41) of the space of the representation. 

Twisted representations: The map 11^-2 cannot be realized by any similarity transformation 
on the space of the representation. 

Indeed, the distinction between these two cases can be simply understood by the identities: 

r „(U„,) = (QnUn^n) /ul, F „(v„) = (fi„V„fi„) /f^. (C.13) 

which tell us that for untwisted representations the local maps Fn leave unchanged the spectrum 
of the generators of the Weyl algebras while this is not anymore true for the twisted ones. 

In particular, these statements imply that for even chain and untwisted representations the SG 
transfer matrix is similar (up to the inhomogeneity flipping) to the r2-transfer matrix while this 
is not anymore the case for twisted representations. 

Remark 6. It is worth remarking that while the spectral problem of the SG model and the 
r2-model are in general different the method here introduced to analyzed the SG spectrum can 
be extended to r2-model. Such an extension is of sure interest as it should allow to get also 
for the r2-model the strong statements of the simplicity and completeness of the transfer matrix 
spectrum. 
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